Baryon-to-Dark Matter Ratio from Random Angular Fields 



CO 

o 

CN 

C 

c3 



43 

i , 

Or 
D ■ 
43 ' 



CN 
> 
CN 



CN 



X 



John McDonalcfl 

Lancaster-Manchester-Sheffield Consortium for Fundamental Physics, 
Cosmology and Astroparticle Physics Group, Dept. of Physics, University of Lancaster, Lancaster LAI 4YB, UK 

We consider the baryon-to-dark matter ratio in models where the dark matter and baryon densities depend on 

angular fields 9j and 6^ according to p t j 8™ and p/, oc 8?, with all values of 6^ and 6^, being equally probable 
in a given randomly-selected domain. Under the assumption that anthropic selection depends primarily on the 
baryon density in galaxies at spherical collapse, we show that the probability density function for the baryon- 
to-dark matter ratio r = Q.b /SIdm ls purely statistical in nature and is independent of anthropic selection. We 
compute the probability density function for r as a function of a and p and show that the observed value of the 
baryon-to-dark matter ratio, r « 1/5, is natural in this framework. 

PACS numbers: 



I. INTRODUCTION 

The ratio of the baryon and dark matter mass densities is 
observed to be Q.b/Qdm ~ 1/5, where DM denotes cold dark 
matter [ 1 ] . A number of different strategies to have been pur- 
sued to account for this. A common origin for dark mat- 
ter and baryons, with a conserved quantum number directly 
linking their number densities, leads to asymmetric dark mat- 
ter with dark matter particle mass is in the range 1-10 GeV. 
Alternatively, the dark matter and baryon density could be 
determined by physically similar processes. For example, a 
weak-strength annihilation process might separately account 
for both the baryon asymmetry and thermal relic WIMP dark 
matter densities ||2|-|4j] . 

All of these approaches account for the baryon-to-dark mat- 
ter ratio via a physical mechanism which directly determines 
the ratio as a function of particle physics model parameters. 
The alternative is anthropic selection. This relies on the exis- 
tence of superhorizon-sized domains with varying Q.b/Qdm- 
One way this can be achieved is if the baryon or dark matter 
densities are determined by a scalar field which is massless 
during inflation. An example based on axion dark matter was 
suggested in [5], while a varying baryon density and the as- 
sociated baryon isocurvature constraints were explored gen- 
erally in |0], with a particular example being Affleck-Dine 
baryogenesis. The possibility of superhorizon-sized domains 
of varying baryon number was first proposed in 0, in the 
context of Affleck-Dine baryogenesis, where it was also noted 
that a range of baryon density at spherical collapse was likely 
to be anthropically preferred. 

However, simply varying one or other of the densities can- 
not by itself completely explain why Gig is similar to SIdm- 
For example, in |5|] the dark matter density varies but the 
baryon density is assumed to be equal to that in the observed 
Universe. The argument of focuses on how increasing the 
dark matter density will affect the baryon and dark matter den- 
sities in galaxies, resulting in conditions which are likely to 
disfavour the existence of observers. This can explain why 



Q.b/Qdm ~ 1 only if we can understand why the observed 
baryon density is close to some critical density above which 
the evolution of observers is disfavoured. But how did the ob- 
served baryon density in our domain get close to this critical 
density in the first place? To explain this we need to let both 
the baryon and dark matter densities vary between domains. 

Here we consider a plausible class of models for the gener- 
ation of domains of varying baryon and dark matter densities. 
We will assume that both the baryon and dark matter den- 
sities are determined by separate angular variables, derived 
from fields which are effectively massless during inflation. A 
specific example is a combination of axion dark matter |i and 
Affleck-Dine baryogenesis dlH. In the following we will 
study this class of model completely generally. The anthropic 
selection criterion we will apply is that there is an anthrop- 
ically preferred baryon density in galaxies at spherical col- 
lapse, such that the probability of the evolution of observers 
is insensitive to the baryon-to-dark matter ratio. This is con- 
sistent with the idea that baryon densities in galaxies which 
are significantly larger than that in the observed Universe are 
anthropically disfavoured di 0] • 

The paper is organized as follows. In Section 2 we intro- 
duce the model and establish the baryon and dark matter den- 
sities in galaxies at spherical collapse. In Section 3 we derive 
the probability density function for the baryon-to-dark matter 
ratio in domains which have equal baryon density at spherical 
collapse. We then apply this function to the case of axion dark 
matter combined with a baryon asymmetry determined by a 
random CP phase (Anthropic Baryogenesis |6[]) and to other 
plausible models. In Section 4 we discuss the plausibility of 
the anthropic selection criterion we have used. In Section 5 
we summarize our conclusions. 



II. DARK MATTER AND BARYON DENSITIES DUE TO 
RANDOM ANGULAR FIELDS 



The class of model we consider is based on the assumption 
that the dark matter and baryon densities are given by 
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where Qd and % are angular fields which are assumed to have 
equal probablility for all values of the angles. The best known 
example of this is axion dark matter. In realistic models the 
behaviour will change as the angles approach 7t, but we will 
consider the cases of interest to have 9,;, 9t < 1. This is ex- 
pected in anthropic selection models, since the anthropically 
selected domains usually have atypically small values of the 
random angles |5l|6[]. 

Due to either random initial conditions at the onset of in- 
flation or quantum fluctuations of the angular fields, super- 
horizon domains of essentially constant 6 ( / and 0& are natu- 
rally generated q. Therefore a horizon-sized Universe at the 
present time will have a random value of (0^,8^), with all 
such pairs of values having equal probability. 

We next consider the growth of perturbations in domains 
with different baryon and dark matter densities. We assume 
that only parameters which vary between the domains are the 
baryon and dark matter densities. We will keep all other pa- 
rameters constant, in particular the magnitude of the primor- 
dial density perturbation. 

In the following it will be important to distinguish between 
mean baryon and dark matter densities in a given domain at a 
given temperature (denoted by pb and pdm) and the mean den- 
sities is galaxy-sized perturbations at spherical collapse (de- 
noted pb s and pdms)- We first consider the effect on galaxies 
of varying (6 ( /,8f,). Varying (Q c /,Qb) will change the mean 
dark matter and baryon densities in the domain at a given 
T. Increasing the mean density will lead to a greatly en- 
hanced density in galaxy-sized perturbations at spherical col- 
lapse. This is because matter-radiation equality then occurs 
at a higher temperature. As a result, significant perturba- 
tion growth (growth proportional to the scale factor), which 
begins approximately at matter-radiation equality, will begin 
at a higher temperature and so the perturbation will become 
non-linear and undergo spherical collapse at a correspond- 
ingly higher temperature and density, enhancing the baryon 
and dark matter densities in galaxies. 

We next derive the enhancement factors for the baryon and 
dark matter densities in galaxies at spherical collapse. Let the 
mean baryon and dark matter densities at a given T in a refer- 
ence domain be p;, „ and p ( / m „. (This could be our domain, but 
we will keep it general.) Suppose in a different domain these 
are enhanced by factors ft, and fd, such that p^ = fbPbo and 
9dm = fdPdm o- The total density at T in this domain is then 



Therefore 



PTOT — fdPdm o + fbPb t 



Ptot _ {ft + faro) 



Ptot, 
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where r is the baryon-to-dark matter ratio in the reference 
domain and we have defined y. The matter-radiation equality 
temperature is determined by 



where K = (% 2 g(T) /9QM 2 p ) l l 2 and g(T) is the effective num- 
ber of massless degrees of freedom. Thus 
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Assuming that the galaxy-sized perturbations enter the hori- 
zon during radiation-domination (this will generally be satis- 
fied except at very small values of r), the temperature at which 
spherical collapse occurs will also be enhanced by y, since the 
rate of growth of the perturbations is the same once the uni- 
verse becomes matter-dominated. The baryon and dark matter 
densities at spherical collapse will then be enhanced by a fac- 
tor y 3 from the earlier time of spherical collapse and by factors 
fd and ft, from the enhancement of the dark matter and baryon 
density at a fixed temperature. Therefore 
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where pt, s and pdms are the baryon and dark matter densities 
in galaxies at spherical collapse 1 . 

The anthropic selection criterion we will apply in the fol- 
lowing is that the probability of the evolution of observers is 
dependent only on the baryon density at spherical collapse and 
is independent of the dark matter density in galaxies and so of 
r. This is not likely to be exactly true, but we will argue later 
that it should be a good approximation to a realistic anthropic 
selection criterion. We will see that the effect of this condition 
is that the probability of r factors out of the total probability 
and becomes purely statistical in nature, allowing us to deter- 
mine its probability density exactly. 

For a given value of pi, s , the value of fd is constrained to 
be a function of ft,, 
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This can also be written in terms of 8^ and 9/,, using / ( / = 
(Qd/®do) a and ft, = {Bb/®bo)^- For simplicity, in the follow- 
ing we will absorb 8,./ and 8/, into the definition of Qd and 



KT eq - PTOT(Teq) , 



(4) 



These expressions are equivalent to the more complicated expression given 
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d b , so that f d = 6" and f b = 9^, with B do =d bo = l. Then 
the constraint is 



and the constraint Eq. ( fTOb to obtain 9/, and 8^ as functions of 
r and p bs 
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The physical interpretation of this constraint is that there are a 
range of different mean baryon and dark matter densities (and 
so baryon-to-dark matter ratios) which can lead to the same 
baryon density at spherical collapse. For example, reduc- 
ing the mean dark matter density causes the matter-radiation 
equality temperature to decrease, reducing the baryon density 
at spherical collapse. This can be compensated by increasing 
the mean baryon density. 
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where we have defined A = (1 + r )/p bSo and we have kept 
the dependence on r and p/„ explicit. The derivatives as a 
function of r and p b s are then 



m. THE PROBABILITY DENSITY FUNCTION FOR THE 
BARYON-TO-DARK MATTER RATIO 

We next calculate the probability density function for r in 
domains satisfying the constraint that the baryon density at 
spherical collapse is between Pb So and p bSo +Api, So , where 
Apbs„ Pbs„ and Api, s is a small change over which the an- 
thropic probability of p b s does not vary significantly. 

Since the probability of each value of 9^ and each value of 
C / is equal, the probability of finding the pair (Qb,®d) in a 
given domain to be within a range of values corresponding to 
an area of the (Q b ,Q d ) plane is simply proportional to the area. 
In Figure 1 we plot the constraint Eq. (TTOb for p bs equal to 
Pbs„ and for p bs equal to p iS[) + Ap/„ o , for the case (a,P) = 
(2. 1). The total probability for p bs to be in the range p bSg 
to Pb s + Ap b s is men proportional to the area between these 
curves. 

To obtain the probability density with respect to r, we need 
to find the area between the curves for a small change Ar. In 
Figure 2 we show the change in (8i,,9 ( /) from a point on the 
first curve due to (i) a change in r for fixed p/, v and (ii) a 
change in p^ s for fixed r. The corresponding area is then the 
area of the parallelogram with sides described by the vectors 
vi and V2, 



Area = |vi xv 2 | = \AQ d 8Q b - AG b dQ d \ 
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To evaluate the derivatives, we impose the relation between 
Q c i, Qb and r, 
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The area is then 
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where 
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Although the first curve is defined to have value p bSo , we have 
now dropped the subscript and consider the probability den- 
sity to be a function of Pb s - The terms K and Api ?s are inde- 
pendent of r and so can be absorbed into a normalization for 
the probability density as a function of r. Thus the probability 
density for a domain with a fixed value of pf, s to have a given 
value of the baryon-to-dark matter ratio r is 



f{r)=N\- 
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where N is a normalization constant which depends on a and 
P- 

The integral of Eq. (1231 from r = to °° is convergent pro- 
vided that a > (convergent as the upper limit tends to °°) 
and a > (3/3 (convergent as the lower limit tends to zero). 
These conditions are likely to be satisfied for realistic models, 
in which case the probability density is well-defined. 

From Eq. (ETT i we see that the probability density for p^ s 
and r factors into an r dependent factor and apj,, dependent 
factor. The total probability density of finding a given value 
of pb s and r is then obtained by multiplying the probability 
density in Eq. ( f2TT > by an anthropic selection factor which de- 
pends only on pt, s - 

The interpretation of the probability density Eq. (l23~t is that 
it gives the probability for any observer with a given value of 
Pi, s to find themselves in a domain with a given value of r. 
The probability of r is purely statistical, since r is assumed to 
have no anthropic effect on the evolution of observers with a 
given value of pi, s , and simply represents the relative number 
of domains with p/„ which have a given value of r, all of 
which are equally likely for the observer to have evolved in. 
This is justified in so far as the dark matter density in galaxies 
does not have a strong effect on the evolution of observers 
i.e. the baryon density at spherical collapse is the dominant 
anthropic parameter. 



Nevertheless, we find that 65% of domains have r in the 
range 0.1-10. Therefore models of this type are clearly able 
to explain why the observed density of baryons is within an 
order of magnitude of the observed density of dark matter. 

The probability density favours lower values of r within a 
given range. This can be seen in Table 2, where we show the 
range r = 0.1 — 1.1 divided into intervals with Ar = 0.2. This 
shows that domains with r around 0.2 are favoured out of the 
subset of domains with r = 0.1 — 1.1 . 




FIG. 1: Region of the (9^, 6j) plane satisfying the constraint of equal 
baryon density at spherical collapse 



A. Application to the axion-anthropic baryogenesis model 



As a first example, we will apply the probability density 
Eq. (1231 to the case of axion dark matter combined with an an- 
thropic baryogenesis model, for example Affleck-Dine baryo- 
genesis la 01 ■ In this case p^ °< Q a and p^ °e fo , corresponding 
to (oc,P) = (2, 1). The probability density is then 
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N 



r 3/8 (1+r) 9/8 



where the normalization factor is given by 
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In Figure 3 we plot the probability density per logarithmic 
interval, r x /(r), which gives an indication of the likelihood 
of finding r close to a given value. (f(r) itself is divergent as 
r — > 0.) The maximum is at r = 1 .25. 

In Figure 4 we show the probability of r being within ranges 
which differ by factors of 10. The corresponding probabilities 
are given in Table 1. We find that 30.3% of domains have r in 
the range 0.1-1. Thus it is natural in this model for the baryon- 
to-dark matter ratio to have the observed value r — 0.2. How- 
ever, it is not the case that values similar to 1 are inevitable. 
From Figure 4 and Table 1 we see that 13.5% of domains have 
r < 0.1, with 3.2% having r < 0.01, while 21.8% of domains 
have r > 10, with 9.5% having r > 50. Thus there is also 
a significant probability that a given observer will evolve in 
a galaxy which is strongly dominated by dark matter or in a 
galaxy with almost no dark matter. 
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FIG. 2: Calculation of probability density function. 



B. Other models 

We show in Table 3 the probabilities for three other plausi- 
ble values of a and p: (i) (oc,p) = (1,1), (ii) (oc,p) = (1,2) 
and (iii) (a, p) = (2,2). We see that all of these models can 
naturally account for the observed value of r, with the frac- 
tion of domains in the range r = 0.1 — 1 being 40.9%, 24.9% 
and 47.6% for models (i), (ii) and (iii) respectively. Therefore 
the ability of this class of model to account for the observed 
baryon-to-dark ratio is not strongly dependent upon the values 
of a and p. Model (ii) has a much larger fraction of domains 
with a large dark matter content, corresponding to r < 0.1. 
This can be understood as due to the stronger dependence 
of Qd on p ( [ (B c i oc p d ) and weaker dependence of 9^ on p^ 
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FIG. 3: Probability density function per logarithmic interval as a 
function of r for the case (a, (3) = (2, 1). 
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FIG. 4: Probability of the baryon-to-dark matter ratio for the axion- 
anthropic baryogenesis model, (a, (3) = (2, 1). 



(06 p b ). This means that a greater increase in the prob- 
ability of a domain is obtained by increasing p c / by a factor 
than by increasing Pf, by the same factor. 



r 


P 


0-0.1 


0.135 


0.1-1 


0.303 


1-10 


0.344 


10-100 


0.154 


100 -°o 


0.064 


0-0.01 


0.032 


10-50 


0.123 


50-oo 


0.095 



TABLE I: Probability of r for the axion-anthropic baryogenesis model (a = 
2,0 = 1). 



r 


P 


0.1-0.3 


0.115 


0.3-0.5 


0.073 


0.5-0.7 


0.054 


0.7-0.9 


0.042 


0.9-1.1 


0.035 



TABLE II: Probability of r for the axion-anthropic baryogenesis model (a = 
2, (5= 1) in the range 0.1-1.1. 



r 




P(l,2) 


P(2,2) 


0-0.1 


0.296 


0.633 


0.385 


0.1-1 


0.409 


0.249 


0.476 


1-10 


0.248 


0.101 


0.214 


10-100 


0.042 


0.015 


0.086 


100 -oo 


0.005 


0.002 


0.035 


0-0.01 


0.093 


0.380 


0.189 


10-50 


0.037 


0.013 


0.070 


50-oo 


0.010 


0.003 


0.052 



TABLE III: Probability of r for alternative models with (a, (5) = (1,1), (1,2) 
and (2,2). 



IV. PLAUSIBILITY OF THE CRITERION FOR 
ANTHROPIC SELECTION 

We have assumed that the main parameter determining the 
evolution of observers is the baryon density at spherical col- 
lapse, and that the density of dark matter at spherical collapse 
has only a small effect. How good an approximation is this 
likely to be? 

Suppose we consider the effect of (i) increasing the baryon 
density at spherical collapse relative to our own domain while 
keeping the dark matter density constant and (ii) increasing 
the dark matter density at spherical collapse while keeping 
the baryon density constant. 

We first consider the effect of increasing the baryon density. 
It is believed that the balance between dynamical collapse and 
radiative cooling favours formation of galaxies with baryon 
mass ~ 10 12 M Q , corresponding to the largest mass for which 
a collapsing gas cloud can cool, and that this mass is indepen- 
dent of the conditions at spherical collapse (being expressible 
in terms of fundamental constants) Q8|] . In this case, increasing 
the baryon density will mean that the volume at spherical col- 
lapse which contains ~ 10 12 Mq becomes smaller, therefore 
galaxies will form from smaller length perturbations. As a re- 
sult there will be more galaxies within a given volume than in 
our domain. (After spherical collapse, we can expect the time 
necessary for observers to evolve to be similar in all domains, 
~ 1 — 10 Gyr.) Since stars are forming out of an initially more 
dense gas cloud, it is also plausible that larger and therefore 
shorter-lived stars will tend to form. In addition, since the 
density perturbation is unchanged, these smaller length per- 
turbations, which would correspond to galaxy subhalo scales 
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in our domain, can be expected to be highly clustered. There- 
fore we can expect multiple mergers of galaxies, resulting in 
large densities of stars and potential disruption of planet or- 
bits. Moreover, if the baryon density becomes large enough, 
the size of the volume with ~ 1O 12 M will become smaller 
than the size of observed galaxy discs, in which case the whole 
picture of galaxy formation would have to be radically altered. 
In particular, as noted in |7|], the density of stars will eventu- 
ally be such that the spacing between stars is comparable to 
the size of a solar system, at which point stable planet orbits 
will become unlikely; this would represent the upper limit on 
the baryon density at spherical collapse. Therefore it is plau- 
sible that increasing the baryon density at spherical collapse 
will have a strong effect on the probability of forming long- 
lived stars with stable planet orbits in a radiation-friendly en- 
vironment, which is the primary requirement for the evolution 
of observers. This is consistent with the idea that anthropic 
selection is sensitive to the baryon density at spherical col- 
lapse and that there is an anthropically preferred value of the 
baryon density at spherical collapse. It is also consistent with 
the naive expectation that increasing the baryon density in the 
volume corresponding to our galaxy at spherical collapse will 
result in a correspondingly larger density of stars, which is the 
basis of lH and H. 

We next consider the effect of increasing the dark matter 
density at spherical collapse while keeping the baryon den- 
sity constant. In this case the perturbations corresponding to 
baryon mass ~ 1O 12 M will have the same size as in our do- 
main. The only effect will be to increase the dark matter den- 
sity in the galaxy halos. While this may have some effect 
on the subsequent evolution of the galactic disc, it would not 
be expected to have a catastrophic effect on the probability of 
forming long-lived stars with stable planet orbits, except in the 
limit of a very large increase in the dark matter density, cor- 
responding to r — > 0. (With fixed baryon density at spherical 
collapse, r effectively parameterizes the dark matter density at 
spherical collapse.) However, domains with larger dark mat- 
ter densities, corresponding to r < 0.1, generally account for 
only a small fraction of the total probability for r. For exam- 
ple, in the case of axion dark matter combined with anthropic 
baryogenesis ((a, p) = (2, 1)), we find that only 13.5% of do- 
mains have r < 0.1. We do not expect domains with smaller 
dark matter densities, corresponding to r > 0. 1 , to have a very 
strong anthropic dependence on r. Moreover, if the effect 
of a large dark matter density is to disfavour the evolution 
of observers, the effect would be to the redistribute the to- 
tal probability density as a function of r to larger values of 
r, which only strengthens the conclusion that r = 0.2 can be 
naturally understood in this framework. Therefore our con- 
clusions should be robust with respect to the anthropic effects 
of the dark matter density. 



There have been other approaches to determining the ax- 
ion dark matter density anthropically ll9l-HTIl. In particular, in 
ifioll it was noted that decreasing the temperature at matter- 
radiation equality while keeping the density perturbation con- 
stant implies that galaxies will not be able to achieve spherical 
collapse before entering the dark energy dominated era. This 
imposes a lower bound on the axion density. We have not 
considered dark energy here, but we note that if we are to un- 
derstand the coincidence of the dark energy density with the 
observed baryon and dark matter densities, it is likely that the 
dark energy density must also be domain-dependent, in which 
case anthropic constraints based on a fixed dark energy den- 
sity do not apply. We will return to this issue in future work. 



V. CONCLUSIONS 

We have shown that baryon and dark matter densities which 
depend on random angular fields, when combined with the as- 
sumption that the dominant anthropic selection parameter is 
the baryon density in galaxies at spherical collapse, can eas- 
ily account for why the observed baryon-to-dark matter ratio 
is approximately equal to 1/5. The total probability for the 
baryon and dark matter densities in a given domain is factor- 
ized into an anthropic factor, which depends only on pt, s , and 
a purely statistical factor, which depends only on r. 

The model is crucially dependent upon both the baryon and 
dark matter densities being able to randomly vary between 
superhorizon-sized domains. To make this clear, we can con- 
trast our model with the model of JD], in which only the axion 
density varies, with the baryon density being fixed to its ob- 
served value. In this case changing the dark matter density 
changes the baryon density at spherical collapse. Therefore 
r is determined anthropically, not statistically as in the model 
we have presented here. Moreover, in the model of [HI, the 
baryon density in our domain must be tuned to be close to 
some anthropically critical baryon density, so that increases 
in the dark matter density to values larger than we observe are 
anthropically disfavoured. No such tuning is required in the 
model presented here, since the baryon density can be natu- 
rally anthropically selected, and the value of r is determined 
purely statistically and not anthropically. 
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